Dirichlet to Neumann Maps 
oi for Infinite Quantum Graphs 

(N 

Robert Carlson 
Department of Mathematics 
University of Colorado at Colorado Springs 
rcarlson@uccs.edu 

<3 ■ September 15, 2011 

43 ' 

Abstract 
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■ for elliptic equations on a large collection of infinite quantum graphs. 
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1 Introduction 



The recent surge of activity in analysis on graphs is wide ranging, encom- 
passing the spectral theory of finite graphs [TTJ [13], physics inspired prob- 
lems on finite or infinite quantum graphs [32] , and resistance network models 
pH EDI [22], often related to probability [15 ] 129]. to mention a few high- 
lights. The study of harmonic functions is a common theme, particularly 
when the work is related to probability The demands of harmonic func- 
tion theory and probability have inspired studies of infinite tree boundaries 
[T0[ IT2] and some work on boundaries for more general graphs [TBI [20], [30] . 

Extending the previous work beyond harmonic functions and trees, this 
paper treats elliptic boundary problems, in particular the Dirichlet to Neu- 
mann map, for a large class of infinite quantum graphs. Broadly speaking, 
Dirichlet to Neumann maps describe the relationship between the value of a 
function / : dB — > R on the boundary OB of some spatial domain B and the 
normal derivative at the boundary of an extension u : B — > R of /. Typi- 
cally, the function u satisfies a differential equation in B, with solutions being 
uniquely determined by the boundary values /. Physically, this formulation 
is used to describe the current flowing out of a domain in response to an 
applied boundary voltage, or the heat flux at the boundary in response to a 
fixed boundary temperature distribution. 

Calderon [5] inspired substantial work on the inverse problem of deter- 
mining (nonconstant) interior electrical conductivities in a domain in R^ 
from the voltage to current map at the boundary; see [28] for an overview. 
Related Steklov eigenvalue and expansion problems are treated in [TJ. The 
Dirichlet to Neumann map and related problems have also been studied for 
a variety of network models. Classical resistor networks are considered in [PJ. 
A quantum graph Dirichlet to Neumann map was used in [I] to develop a 
finite tree version of the Borg-Levinson inverse eigenvalue theorem. Bound- 
ary control methods as developed in [2] [3] provide an alternative technology 
for Dirichlet to Neumann mapping problems. 

The class of quantum graphs treated here is motivated by the challenge 
of modeling enormously complex biological networks such as the circulatory 
system [6] , the nervous system [26] , or the pulmonary network [251 EI] • The 
graphs Q are typically infinite, but have finite diameter and compact metric 
completions Q. They also satisfy a 'weakly connected' condition introduced in 
[7] . Not only do these conditions provide productive idealizations of biological 
networks, but they also appear in the study of continuous time Markov chains 
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with 'explosions' [8]. 

This work resolves a number of distinct problems while extending the 
earlier work on harmonic functions on infinite trees to encompass elliptic 
boundary problems for a large class of infinite quantum graphs. The results 
include an existence theorem for the Dirichlet problem for a class of ellip- 
tic equations. The subsequent study of the Dirichlet to Neumann map is 
motivated by the problem of discussing flows through the graph boundary. 
Classical work on the Dirichlet to Neumann map is usually handled in the 
context of spatial domains with smooth boundaries. Part of the challenge 
here is that the boundaries dQ of infinite graph completions usually lack any 
different iable structure, so alternative formulations are required. 

The second section begins with a review of developments from [TJ. Prop- 
erties of weakly connected metric graph completions and their algebra A of 
test functions are recalled. These graph completions have totally discon- 
nected boundaries, so there is a rich collection of clopen sets, that is sets 
which are both open and closed, which play an important role in the analy- 
sis. The elliptic equations u" = qu with q > and standard quantum graph 
vertex conditions are then introduced. A maximum principle is established 
for these equations; it proves to be a key tool in the analysis. A solution of 
the Dirichlet problem for harmonic functions developed in [7j is extended to 
include the more general elliptic equations and a larger class of graphs. 

The third section introduces the Dirichlet to Neumann map, beginning 
with finite graphs, where classical derivatives at the boundary are available. 
Computations for the a — (3 tree illustrate the extraordinarily complex behav- 
ior of derivatives of simple harmonic functions on infinite graphs. Using the 
voltage to current map physical interpretation for guidance, it proves pro- 
ductive to think in terms of current flow through a boundary set rather than 
current flow at a point. This idea is first implemented using test functions 
from A to define a Dirichlet to Neumann function A g (F, Q). Here F may be 
any continuous function on the boundary dQ of the graph completion, but 
Q is restricted to be a clopen subsets of dQ. 

To ensure more regularity of the Dirichlet to Neumann function with 
respect to the set Q, the fourth section considers the case when F is the 
characteristic function of a clopen subset of dQ. The restricted functions 
A q (F, fl) extend in a standard way to signed Borel measures on dQ with the 
expected positivity when F = ln(i) and Q C Q(l). In the final section, an 
extension by linearity in the first argument yields a Dirichlet to Neumann 
map which is a densely defined nonnnegative symmetric operator from the 
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continuous functions on dQ to the dual space of Radon measures. 

2 Foundations 

This section reviews some of the terminology and results from [7J, which 
should be consulted for proofs and additional information. After introduc- 
ing the geodesic metric, weakly connected metric graph completions are de- 
scribed. These graph completions have totally disconnected boundaries, with 
a rich collection of clopen sets, which are both open and closed. Weakly con- 
nected graph completions have a closely associated algebra A of 'eventually 
flat' functions. Following this review, the elliptic equations u" = qu with 
q > and standard quantum graph vertex conditions are introduced, and 
a maximum principle established. A solution of the Dirichlet problem for 
harmonic functions developed in [7] is extended to include the more general 
elliptic equations and a larger class of graphs. 

2.1 Metric graphs 

Q will denote a locally finite graph with a countable vertex set. Interior 
vertices are those with more than one incident edge, while boundary vertices 
have a single incident edge. Q is assumed to be connected, but connectivity 
is not assumed for subgraphs that appear in the course of proofs. 

To study differential equations on Q, edges e n in the edge set £ are iden- 
tified with real intervals of finite length l n . Loops and multiple edges with 
the same vertices are not directly considered, but can be accomodated by 
adding vertices. In the usual quantum graph style, the Lebesgue measure 
for intervals may be extended to a measure for Q. The Euclidean length on 
the intervals is extended to paths of finitely many intervals by addition, and 
then to a 'geodesic' distance between points p\ and P2 by taking the infimum 
of the lengths of paths joining p± and pi- 

As a metric space, Q has a completion Q which is assumed to be compact 
in this work. Without much difficulty [3 Prop. 2.1] one sees that Q is 
compact if and only if for every e > there is a finite set of edges S = 
{ei, . . . , e n } such that for every y £ Q there is a edge G S and a point 
Xk € tk such that d(xk, y) < e. Define the boundary dQ of Q to be the set of 
boundary vertices of Q . The interior of Q or Q is Q \ dQ. The boundary dQ 
of the completion Q is defined to be the union dQ = dQ U [Q \Q]. Since dQ 
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is closed, it too is compact. 

Graphs and their completions are more amenable to analysis when they 
satisfy an additional condition. Say that Q is weakly connected if for every pair 
of distinct points x,y G Q, there is a finite set of points W = {wi, . . . , wk} C 
Q separating x from y. That is, there are disjoint open sets U, V with Q\W — 
U U V, and with x & U , y G V . An alternative characterization is that there 
is a finite set of edges Wg from Q such that every path from x to y contains 
an edge from Wg. 

Completions of trees are weakly connected. Less obviously, so are the 
completions of graphs with finite volume [7J. The following result, similar to 
one in [5], shows that trees can be modified substantially without losing a 
weakly connected completion. 

Theorem 2.1. Suppose the graph Q is obtained from a tree T by adding a 
sequence of edges e n whose lengths l n satisfy lim n / n = 0. Assume there is a 
positive constant C such that 

dg(x, y) < d r (x, y) < Cd g (x, y), x,y G T. 

Then Q is weakly connected. 

Proof. We haven't changed the set of Cauchy sequences of vertices, so Q\Q = 
T\T. Suppose x and y are distinct points in Q \ Q. Let E be an edge with 
length L on the path from x to y in T. Let U x be the connected component 
of T\ E containing x, and let V y be T\ (U x U E). Find iV so that l n < L/C 
for n > N, and for n < N remove edges e n , along with E, from Q, leaving 
Q\. Suppose there is a path in Q\ connecting U x to V y . Then there must 
be vertices u and v separated by E in T but joined by an edge e n whose 
length, by assumption, is smaller than L/C. Then we have dr{u, v) > L, but 
dc(u,v) < L/C, a contradiction. Thus cutting E together with ej, . . . , ejy 
provides a separation of x and y. 

□ 

If Q is weakly connected, then dQ is totally disconnected, meaning that 
connected components of dQ are points. Totally disconnected compact metric 
spaces have a rich collection of clopen sets, which are both closed and open. 
In fact [121 p. 97], if Q is a totally disconnected compact metric space, then 
for any x G Q and e > 0, there is a clopen set U which contains x and is 
contained in the e ball centered at x. 
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The connectivity properties of Q are closely related to the algebra A 
of 'eventually flat' functions (j) : Q — > R which are continuous, infinitely 
differentiable on the open edges of Q, and for which <fi'(x) = for all x in 
the complement of a finite set of edges, and in a neighborhood of each vertex 
v G Q. With pointwise addition and multiplication, A is an algebra. Since 
the constant functions are in A, the following result Lemma 3.5] may be 
combined with the Stone- Weierstrass Theorem to show that if Q is weakly 
connected and compact, then A is uniformly dense in the real continuous 
functions on Q. 

Lemma 2.2. Suppose Q is weakly connected. If Q and Q± are disjoint com- 
pact subsets of Q, then there is a function <f> e A such that < <f> < 1, 

^(x) = 1, x e O, (j)(y) =0, ye 
2.2 Forms and operators 

Suppose q : Q — > R is nonnegative and locally integrable. The main focus of 
this work is the study of certain solutions of the equation 

u" = qu, u-.g^R. 

The equation is to hold in the sense that v! is absolutely continuous, and the 
equation holds a.e. on the edges of Q. In addition the function u is assumed to 
satisfy the standard continuity and derivative junction conditions at vertices 
v with degree at least 2 in Q, 

lim Ui(x) = lim Uj(x), e(i), e(j) ~ v. (2.1) 

rrGe(i)— >v x£e(jr')— >v 

^2 d v Ui(v) = 0. 

Here the derivative d u Ui(v) = u' is computed in local coordinates identifying 
[a, b] with the edge e{i) of length b — a, and a is the coordinate value for v. 
Functions u : Q — > R satisfying the equation and vertex conditions will be 
called g-harmonic functions. 

Given an edge weight function oj : £ — > (0, oo) constant on each edge, 
define the weighted inner product 

(f,g) u = ^2 f( x )9{x)u{e) dx, 

e Je 
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and corresponding Hilbert space L^- Let V^ tU = AnL"^. Introduce the form 
Q(f,9) = I f'(x)9'{x) + q(x)f{x)g(x), (2.2) 



jg 

with domain T>a,w, and the differential expression 

Kf=-^r[-f(x)+q(x)f(x)]. 
w(e) 

For f,g E T>a jUI , integration by parts leads to the alternate form 

W,2) = E / -^-A-f"(x) + q(x)f(x)}g(x)oo(e) = J2 f[L u f(x))g(x) U (e). 

These computations are summarized in the following proposition. 
Proposition 2.3. The bilinear form 

Q(f,g), f,g£i>A,uj, 

is densely defined in L^, has a nonnegative quadratic form Q(f, f), and sat- 
isfies 

Q(f,g) = (Luf,g)u = (f,L u g) u . 

The form Q has an associated inner product [2H pp. 308-318] 

(f,g)u,i = (f,g)u + Q(f,g), 

with the corresponding norm ||/|| U) i. Completing the domain with 
respect to this inner product, construct the Hilbert space 

Integration from x to y over a path 7 of length at most 4d(x, y) yields 
the estimate 

\f(y)-f(x)\ 2 = \ [ f'(t)dt\ 2 < I ldt I \f'(t)\ 2 dt<4d(x,y)Q(f,f). (2.3) 

Since Q is connected with finite diameter, the arguments of [7] show that 
functions / in H 1 ^) are continuous on Q. There is a constant C such that 

sup |/(a;)| <C\\f\\ u , u 
x&g 



and / satisies a Lipshitz estimate 

\f(y)-f(x)\<2d(x,y) 1 / 2 \\f\U 1 . (2.4) 

Consequently, / has a unique continuous extension to Q. 

Let Vq denote the set of functions / e A whose support is contained in 
the union of a finite collection of edges, and such that f(x) = for all x in a 
neighborhood of any boundary vertex. Let Sq denote the symmetric nonneg- 
ative operator acting by on the domain V in L^. So has a self adjoint 
Friedrich's extension £ acting on Lj. If q is not bounded, the description of 
the domain of So is more delicate than described. The reader should either 
consult [2TJ PP- 343-346] for the necessary modifications, or simply assume 
that q is bounded. 

Suppose u)(Q) = X]e w ( e ) < 00 • If / £ A) with ||/|| w = 1, then there is 
some point y G Q with \f{y)\ 2 > u;~ 1 (Cr)- If the connected graph Q has a 
boundary vertex or has infinitely many edges, then there will be some x E Q 
with f(x) = 0. Then Q gives 

u;- 1 ^) < \f(y)\ 2 = \f(y) - f(x)\ 2 < ARQ{f, /), 

leading to the next result. 

Proposition 2.4. Suppose oo(Q) = *Yl e u{e) < oo and dQ ^ 0. Then the 
self adjoint operator Cq on has a strictly positive lower bound. That is, 
there is a C > such that for all f in the domain of C with \\f\\ u = 1, 

2.3 g-harmonic functions 

There is a maximum principle for g-harmonic functions. 

Lemma 2.5. Suppose q is a nonnegative locally integrable function, and 
u : Q — )■ K is q - harmonic. If u has a positive global maximum at an interior 
point of Q, then u is constant. 

Proof. Arguing by contradiction, suppose u is nonconstant with a positive 
global maximum at the interior point £2 £ Q- Find another interior point 
Xq with u(xq) < u{x2). The first step is to find an interior point X\ with 
u{x\) = u(x 2 ) such that every neighborhood of x\ contains points x with 
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u(x) < u{x\). Recalling that Q is connected, let 7 : [0, 1] — > Q be a path in 
the interior of Q with 7(0) = Xq and 7(1) = Xi- If 

to = sup {t G [0, 1] I ^(7(5)) < u(x2) for all s < t}. 

then xi = 7(^0) has the desired properties. 

If Xi is not a vertex, then u'(xi) = by Calculus. If Xi is an interior ver- 
tex, then u' n (xi) = for each edge e n incident on xi by the vertex conditions 

The function u" is locally integrable, so [27J P- HO] for x near xi, 

PX PX 

u'(x) = / u"(i) dt = q{t)u(t) dt > 0. 

J XI J Xl 

Since u is not constant in any neighborhood of xi, u'(x) > for some x 
arbitrarily close to x\. Since u' is continuous on each edge of Q, 

PX 

u(x) — u[x}) = / u'(t) dt > 0, 

J Xl 

so Xi is not a local maximum. □ 

Lemma 2.6. With the hypotheses of Lemma \2.5\ suppose u : Q — > K is 
continuous, nonconstant, q-harmonic on Q , and < u(x) < 1 for all x e dQ . 
Then < u(x) <1 for allx eG\ dQ. 

Proof. By Lemma [2.51 u(x) < 1 in Q \ dQ. Similarly, if u(xq) < for some 
Xq G Q, then — u, has a positive interior maximum, which is not possible. 

This leaves the case when u(xo) = for some xo G Q \dQ. Interior points 
of edges may be treated as vertices of degree two, so our focus is on a vertex 
v with degree at least two, with u{y ) = 0. Let Ui denote the restriction of 
u to an edge e« incident on v. Since Ui{x) > 0, the derivatives, computed in 
outward pointing local coordinates, satisfy d u U{(v) > 0. The vertex condi- 
tions (12. ip then force u'^v) = for each incident edge e^. On the function 
Ui satisfies the ordinary differential equation u'[ = qu^, with initial condition 
u{v) = u'(v) = 0, so u(x) = for all x G e^. These observations also imply 
that u is identically along any path once it vanishes at any point in Q\dQ . 

□ 

Extending Theorem 4.2 of [7] with a similar proof, the next result shows 
that if u~ 1 q G L 2 (u), then the Dirichlet problem is solvable for (/-harmonic 
functions. In particular the Dirichlet problem for the usual harmonic func- 
tions (q = 0) is solvable whenever Q is weakly connected and compact. 
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Theorem 2.7. Suppose Qis connected, while Q is weakly connected and com- 
pact. Assume there is a weight u such that u(Q) = J2 e u(e) < oo, and the 
nonnegative locally integrable function q satisfies 

I {uj^qfw < oo. (2.5) 
Jg 

IfdQ 7^ 0, then every continuous function U : dQ — > K. has a unique extension 
to a continuous function u : Q — > K that is q-harmonic on Q . 

Proof. As noted in the remarks preceeding Lemma 12.21 A is uniformly dense 
in the real continuous functions on Q. After extending U to a continuous 
function on Q, find a function g n G A with 

max \ U{x) — g n (x)\ < 1/n. 
xeg 

Since g n G A, the function D 2 g n is bounded with support contained in the 
union of a finite collection of edges. Since g n is bounded, (12. 5 p implies that 

uj- 1 [D 2 -q]g n eL\u). 

The operator £ is invertible, so the function 

K = C q 1 uj- 1 [D 2 -q]g n (2.6) 

satisfies 

u- 1 [-D 2 + q}h n = uj- 1 [D 2 -q}g n . 

Since g n G A and h n is in the domain of Cq, u n = h n + g n satisfies the vertex 
conditions (12. ip . In addition h n (x) = for x G dQ. That is, u n is g-harmonic 
on Q with u n (x) = g n (x) for x G dQ. 
Suppose m <n. Then 

2 

max kt m — M n | < — . 

x&ag m 

Since both ±(u m — u n ) are g-harmonic on Q, the maximum principle implies 

2 

max \ u m — u n \ < — , 

xeg m 

so {u n } is a uniformly Cauchy sequence on Q, with a continuous limit u. 
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On each edge e the solutions of u" = qu form a two dimensional vector 
space. Recall [2H pp. 4-6] that any two norms on a finite dimensional 
vector space induce the same topology, so the uniform convergence to of 
\u m (x) — u n (x)\ for x G e implies the uniform convergence to of \u m (x) — 
u n (x)\ + \u' m (x) — u' n (x)\. Recasting solutions of u" = qu as solutions of 
an integral equation and taking limits shows that u(x) = lim n u n (x) is the 
desired g-harmonic function. 

□ 

3 Introducing the Dirchlet-Neumann map 

Beginning with finite quantum graphs, where classical derivatives at the 
boundary are available, this section introduces the Dirichlet to Neumann 
map. Sample computations for a family of infinite a — (3 trees demonstrate 
the complex behavior possible for derivatives of harmonic functions on infi- 
nite graphs. For continuous boundary functions U and clopen sets Q C dQ, 
a preliminary Dirichlet to Neumann function A q (U, Q) is defined using test 
functions from A. Since the existence of g-harmonic extensions is fundamen- 
tal, the hypotheses of Theorem 12 . 71 are assumed. Moreover, q is subsequently 
assumed to be integrable over Q. 

3.1 Finite Graphs 

In case Q is a finite connected metric graph, dQ is just the set of vertices 
of degree one. The hypotheses of Theorem 12.71 are satisfied if oj{e) = 1 for 
each edge e, so for any function U : dQ — > R there is a unique q - harmonic 
extension u : Q — > KL The function u will have an absolutely continuous 
derivative. For each boundary vertex v, let d u denote the derivative computed 
with respect to local coordinates pointing outward at v. The Dirichlet to 
Neumann map A q acting on the vector space of functions on dQ is defined 
by 

A q U{v) = d v u{v), v e dQ. 

For a finite graph, introduce the Hilbert space l 2 (dQ) of functions F : 
dQ — > R, with the usual inner product 

(F,G) B0 = y £ i F(y)G(v). 
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Extending functions F : dQ — > M to smooth functions / : Q — > K. satisfying 
the standard interior vertex conditions (I2.ip . integration by parts gives 

/ ("/" + qf)g = - E (&/ («) + / f'9 + Qf9- (3- 1 ) 
^ e vedg ^ g 

With L/ = — /" + qf, a second integration gives 

f (Lf)g = lf(v)d»g(v) - g{v)d v f{v)\ + / f(Lg). (3.2) 

If L/ = and g = 1, (13. 2p reduces to 

E = / /(^)- = / ( 3 - 3 ) 

When q = this equation corresponds to 'conservation of current'. In case 
Lf = Lg = 0, (13. 2p reduces to 

(A q F,G) 9g = 9(v)d v f\v) = /(«)M«) = (^A ? G) a g. 

If, in addition, G = F, then (13. ip gives 

(A q F,F) 9g = [ (f'f + qf. (3.4) 

Finally, if A q F = 0, then (13.41) with q > shows that /' = everywhere, 
which implies / is constant on Q, and q = almost everywhere. These basic 
computations are summarized in the following proposition. 

Proposition 3.1. For a finite graph Q with nonempty boundary dQ, the 
Dirichlet to Neumann map A q is self adjoint and nonnegative on l 2 (dQ). 
The null space of A q is {0} unless q = almost everywhere, in which case 
the null space consists of the constants. 

Using the g-harmonic extensions / of F and g of G, (13. ip expresses the 
bilinear form (A q F, G)g g as an integral over Q, 



(A q F,G) dg = [ f'g' + qfg. 
Jg 
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It is also possible to express A q F using test functions in (13. 2p . Let <p(x) be 
a smooth function with <fi' e {w) = for all edges e and vertices w G Q. If 
</>(i>) = 1 at the boundary vertex v, while <p(w) = at all other boundary 
vertices, then (13.21) gives 

A q F(v) = d v f(v) = [ fL<j>. (3.5) 
Jg 

3.2 Infinite network example: the a — (5 tree 

Since the derivatives of harmonic functions on infinite graphs can exhibit 
complex behavior, extending the Dirichlet to Neumann map from finite to 
infinite networks will require new ideas. To illustrate this complexity, con- 
sider an infinite tree T a , as illustrated in Figured! T a has a root vertex t> 
of degree 1 with adjacent vertex v\. All vertices except vo have degree 3. 
Organize the vertices into combinatorial depth levels, so v is at level 0, Vi 
is at level 1, and if v n is at level n, its children are at level n + 1. 

vO 



vl 




Figure 1: An a — (3 tree 



The edge weights are described by a scaling parameter a, where a and 
= 1 — a satisfy < a, (3 < 1. For any vertex v G T other than vo, label the 
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single incident parent edge e p (v), and the two child edges e a (v) and ep(v), 
with lengths l p , l a and lp respectively. The lengths satisfy l a = al p , lp = f3l p . 

Given a linear function u(x) = cx + d on the edge [vq, Vi], extend u to T a 
by requiring u to be continuous on T a and linear on each edge, with 

a + p 
— cx 

u \ x ) = — n^ u (*)> * G e p( u )> x G e ^ v )- 

a + p 

The derivatives are computed in local coordinates pointing outward from v.. 
One checks easily that u is harmonic on T a , has the same value at all vertices 
on the same level of T, and has a finite limit as the level index n — > oo. 

Paths descending through the levels are determined by sequences {cr n } 
with a n G {a, (3}. The derivative below the vertex vn at level N will be 

u'(v+) 

Note that limjv^oo u '( v n) = independent of the path, but the decay behav- 
iors are path dependent and diverse. 

3.3 Test function formulation 

As noted in (13.51) . there is a test function formulation of the Dirichlet- 
Neumann map for finite graphs. For infinite graphs the points in dQ are 
generally not isolated, so the previous formulation is not applicable. An 
alternative is available when Q is compact and weakly connected, with func- 
tions in A serving as test functions and clopen subsets of dQ generalizing 
boundary vertices. 

The set dQ is a closed subset of Q, and dQ is totally disconnected if 
Q is weakly connected. Since dQ is a totally disconnected compact metric 
space, every open ball B e (x) contains a clopen neighborhood of x. Apply- 
ing Lemma I2.2[ the next lemma guarantees a plentiful supply suitable test 
functions. The succeeding lemma describes their features. 

Lemma 3.2. Suppose Q is compact and weakly connected. If Q is a clopen 
subset of dQ , then there is a function G A such that < < 1 and 

00) = 1, xeQ, 0(|/) = 0, yedQ\Q. 



a/3 

w'Or) n — ' u '( v i) = °- 



n=l 



a 



II 
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Lemma 3.3. If (f) G A, then the restriction : dQ — > K /ias finite range R. 
The set {x\<fi(x) ^ i?} zs contained in the union of finitely many edges. 



Proof. Let E be a nonempty finite set of closed edges containing all points 
x with <p'{x) 7^ 0. For any point z G dQ \ E there is a path 7 : [0, 1] — > Q 
with 7(0) G £, 7(1) = z, and 7(f) G 6? for < t < 1 by Prop 2.3]. 
Since d(z,E) > there is a smallest number t with < t Q < 1, and such 
that j(t) $l E for all t > to- Since is the union of a finite set of closed 
edges, 7 (to) is an endpoint of one of these edges, and 4>(z) = 0(7(t o )), since 
(//(7(t)) = for all t > to- This shows that the restriction : dQ — > R has 
finite range i?, and that {x\<p(x) ^ R} C E. □ 

Assume the hypotheses of Theorem 12.71 are satisfied. Before treating the 
Dirichlet to Neumann map, a preliminary step is to consider the real valued 
Dirichlet-Neumann function A q (U,Q), defined for continuous functions U : 
dQ — > R and clopen sets Q C 9^. The value of A q (U, ft) plays the role of 
the integral of the normal derivative of the g-harmonic extension u over set 
Q C dQ. Using (13. 5p as a guide, suppose <fi E A satisfies 

<f>(x) = l, xeQ, <f>(x) = 0, xedQ\Q. 

The Dirichlet-Neumann function A q (U, Q) is defined by the formula 

A q {U,0) = [ uL<P, (3.6) 
Jg 

The next lemma shows that the dependence on cf) is illusory. 
Lemma 3.4. Suppose Q C dQ is clopen, and that (f>, ij) G A with 

cj)(x) = ip(x) = 1, x G fi, 0(x) = ^(a;) = 0, 3; 6 95 \ fl. 

Then 



uL<fi = / 
g Jg 

Proof. The integration by parts formula (13.21) will be used for certain finite 
subgraphs Q\ of Q. The validity of this formula is based on the assumption 
that the interior vertex conditions (12. ip hold at vertices of degree greater 
than one. A subgraph Q x may include vertices v with degree di(v) in Qi, but 
with strictly larger degree d(v) in Q. To resolve this problem, the boundary 
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dQi is understood in a relative sense, comprising the set of vertices v & Gi 
with d\{v) = 1 or with diiy) < d{y). In the formula ( I3.2p . such a vertex v is 
treated as d\{v ) boundary vertices of degree one. 

Using Lemma 13.31 and the fact that <ft, ip 6 A, there is a finite connected 
subgraph Go ^ G such that for all connected finite subgraphs Gi with Go C 
Gi C G we have 

o = = if/( x ) = <p(x) - ij>(x), xeg\g x , 

and 

= 4>'(x) = ip'(x) = 4>(x) — ip(x), x e dQ\. 

Since Lu = 0,the formula f)3.2p gives 

/ uLcj) = <j){v)d v u{v) = ijj(v)d v u{v) = / uLip. 

This equality gives 



uL(j> — uLifj | = | / uLcj) — uLip\ = | / uq(4> — 
g Jgxs! Jg\g! 

Since u(4> — ip) is bounded and q is integrable, the right hand side can be 
made arbitrarily small by enlarging Q\. □ 

Having established a Dirichlet-Neumann function A q (U,Q), a couple of 
preliminary observations are in order. The first concerns nontriviality. 

Proposition 3.5. IfU : dQ — > K. is nonnegative, continuous, and somewhere 
positive, and if q(x) is integrable and strictly positive, then A q (U,dQ) > 0. 

Proof. By Lemma [2761 the q harmonic extension u of U is everywhere nonneg- 
ative. For the clopen set Q = dQ simply take <fi to be the constant function 
1. Then 

A g (U,dG)= [ uq>0. 
Jg 



□ 



For the second observation, suppose the hypotheses of Lemma 13.41 hold. 
Assume that {Gn} is an exhaustion of Q by a sequence of finite subgraphs, 
so that 

Gn C Gn+l, G = [J Gn- 
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For all n sufficiently large 



<P'(x) = o, 0(x)e{o,i} 



x g dQ,, 



Using the relative notion of graph boundary, let B n = _1 (l)n<9(? n . Applying 
( 13.21) as in the proof of Lemma [33 leads to an approximation of A q (U, Q) using 
finite graph Dirichlet to Neumann maps, 



4 Characteristic functions 

This section focuses on the Dirichlet to Neumann function A q (U,Qi) when 
U = 1q : dQ — > R is the characteristic function of a nonempty clopen set, 



Let denote the g-harmonic extension of 1^ to Q. By Lemma [231 if Vl ^ dQ, 
then the inequality < unit) < 1 holds for all t G (Q\dQ). 

For guidance, consider a finite graph Q and point x G Q with u(x) = 1. 
The Mean Value Theorem together with the equation u n = qun > imply 
that u' n (x) > 0. For y G (dQ \ ft) a similar argument gives u'(y) < 0. This 
section will establish analogous results for infinite graph completions. 

4.1 Regular values 

To provide technical support for an upcoming argument, the notion of regular 
values of a function is adapted to functions defined on networks. Assume that 
h : Q — > M is continuous, with a continuous derivative on the open edges of 
Q. A point x G Q is a critical point for h if x is a vertex or if h'(x) = 0. 
A number y G K is a critical value for h if h~ l (y) contains a critical point. 
Points in the range of h that are not critical values are regular values. The 
next result is a simple form of Sard's Theorem. 

Lemma 4.1. The set of critical values of h has Lebesgue measure 0. 




(3.7) 
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Proof. The set of vertices is countable, so the set of values h(v) for vertices 
v has measure 0. There are countably many edges, so it suffices to prove the 
result for a single edge e. Suppose the length of e is l e , and e is identified 
with the real interval [0,l e ]. 

Since h! is continuous, for e > the set U e = {x\ — e < h\x) < e} is 
open. U e may be written as the union of countably many pairwise disjoint 
open intervals (a„, (3 n ) C (0, l e ). With z n = (a n + (3)/2 and a n < x < (3 n , 

h(x) - h(z n ) = j h'(t) dt. 

Thus the image of the n-th interval has length at most (/3 n — a n )e, implying 
h(U t ) is contained in a set of measure at most el e . □ 

Lemma 4.2. Suppose Q is compact, and that c is a regular value of h. As- 
sume there is no x G dQ with h(x) = c. Then h~ l (c) is a finite set containing 
no vertices. 

Proof. If /i _1 (c) were an infinite set, then by compactness there would be 
an infinite sequence {x n } C h~ l (c) converging to a point z, with h(z) = 
h(x n ) = c. Moreover z is an interior point of some edge e. By the Mean Value 
Theorem and the continuity of h', h'(z) = 0, contradicting the assumption 
that c is a regular value for h. □ 

Level sets of h can be used to construct subgraphs of the metric space 
Q with added vertices and subdivided edges. Let c±,C2 be regular values 
of h with ci < C2, and assume that h(x) ^ [ci,C2] for all x G dQ. Let 
C = /i _1 (ci) U /i _1 (c 2 ), which is finite by Lemma [4.21 Add vertices of degree 
two at the points of C, and subdivide the corresponding edges to obtain a 
modified graphical description of Q. Let Q(ci, c 2 ), the level set graph, be the 
set of points x G Q such that every path from x to dQ contains a point in C . 

Lemma 4.3. If Q is compact, then Q[c\,C2) is a finite graph whose boundary 
consists of points in C. If v is a boundary vertex of Q(c\,C2) with h(v) = c^, 
then d u h{v) > 0, while if h{v) = c\, then d v h{v) < 0. 

Proof. Since dQ is compact, there is an e > such that h(x) [01,02] if 
d(x,dQ) < e. This implies [7J Prop 2.1] that Q(ci,c 2 ) is contained in the 
union of a finite collection E of edges of the (original) graph Q. Assume that 
E is chosen so that every edge e G E contains some point of Q(ci, c 2 ). 
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Consider the following observations about an edge e G E, identified with 
an interval [a,b]. First, if x ,Xi,x 2 G e with x < x\ < x 2 , and x ,x 2 G C, 
then x\ G {?(ci, C2). Second, if xq, x 2 G e with x < £2, and there is no point 
£1 G C with so < xi < x 2 , then by the Intermediate Value Theorem either xq 
and x 2 are both in G(ci,c 2 ), or neither is. Third, if X\ G C, then h'{x\) 7^ 
implies that an interval containing x\ is a subset of (?(ci, C2). 

Let N be the number of points of C in e. By considering the three cases 
N = 0, N = 1, and N > 2, one finds that the set e fl Q(c\,c 2 ) is a finite 
union of closed edges from the modified graph Q (including points of C as 
vertices). Consequently, G(ci,c 2 ) is a finite graph whose boundary vertices 
are in C. The conclusions about the derivatives at boundary vertices v follow 
from h'(v) ^ 0. □ 



4.2 Boundary measures 

Theorem 4.4. Assume that Q is connected, and that Q is weakly connected 
and compact. Suppose Q,Ei,E 2 are nonempty clopen subsets of dQ with 
Q ^ dQ, E 1 C tt, and E 2 C dQ \ SI. Then 

A q (l n ,E l )>0, A q (l n ,E 2 )<0. (4.1) 

If q(x) = for all x in some neighborhood of dQ , then K q {l^E 2 ) < 0. 

Proof. The two cases in (14. Xp are similar; the E% case will be emphasized. 

By Lemma [221 if «j] is the nonconstant g-harmonic extension of 1^ to Q, 
then 

0<« n (a:)<l, xeQ\8Q. 
Since Q is connected the range of uq is [0,1]. Suppose (f> G A with 

0(s) = l, xeE 1: 00) = 0, xedQ\E 1 . 

The idea of the proof is to use a finite graph C/ > constructed as in 
Lemma 14. 31 to approximate Q. The approximating graph should include 
all points where <fi'{x) 7^ so that 

A g (l n ,E 1 )= / u n L(f>= / u u L(f)+ / u n q4>. (4.2) 

Jg Jg\g 

The integral over ^0 is expressed more transparently using ( 13 . 2 p . The pos- 
sible positivity of the last integral term accounts for the slight difference in 
inequalitites for E\ and E 2 unless q vanishes near dQ. 
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By Lemma 13.31 there is an Si > such that if d(x,dQ \ E{) < Si then 
4>{x) = 0, while if d(x, E{) < Si then <p(x) = 1. The integrability of q implies 
that for any e > there is a S > such that 

/ u n q<P < e, N{S) = {xeG\ d(x, dQ) < S}. 
JN(5) 

Choose S < Si. 

Pick a finite subgraph Qs of Q such that Q \ Qs is a subset of N(S). The 
subgraph Qs may contain boundary vertices of Q, but by adding a finite set 
of new vertices to Q, and subdividing the associated edges, Qs can be chosen 
so that it contains no points of dQ. By Lemma [2T61 the function uq, restricted 
to Qs has a minimum ot\ and a maximum /3i satisfying < a,\ < fii < 1. Now 
pick regular values a and (3 for un so that < ati < a and j5\ < ft < 1. Using 
the previous level set graph selection, define Q = Q(a,fl). Every boundary 
vertex of Qq is in the set N(S). 

There will be a nonempty set of boundary vertices v of Qo with uq(v) = (3 
and <f)(v) = 1. By Lemma [4.31 we then have d u un(v) > 0. Using (13.21) as in 

/ u n L(f) = 9 u un(v) > 0. 

Together with 

< / u n q<p < e, 
JQ\Qo 

the result is established. □ 

Fix a reference clopen set Q as in Theorem 14.41 with the corresponding 
g-harmonic extension : Q — > R of 1q : dQ R. By Theorem 14.41 the 
function A 9 (1q, E) is positive on nonempty clopen subsets E of Q. Consider 
extending the set function A g (ln, E) to a positive measure on fl, and a signed 
measure on dQ. The construction of a measure from a premeasure [TTJ p. 30] 
is used here. Recall that an algebra of subsets of Q is a nonempty collection 
closed under finite unions and complements. The collection of clopen subsets 
of dQ is an algebra. 

Lemma 4.5. Assume that E C £1 is clopen, and {E n } is a countable partition 
of E by clopen sets. Then 

A 9 (l n ,£) = X)A,(ln,£„) 

n 
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Proof. Since Q is compact and E is closed, E is compact. The given partition 
is an open cover of E, so is finite. 

Finite additivity will now follow from (I3.6p . For each set E n there is a 
function (fi n G A such that, for x G dQ, <f> n (x) = 1 if x £ E n and <j> n (x) = 
if x £ E n . The function = ^ n </> n is in A, restricts to the characteristic 
function of E on dQ, and 

A q (l a ,E) = [ un£<f>= [ u n £j2<Pn = J2 A «( 1 ^ E n)- 

□ 

Theorem 4.6. A q (1q,E) extends to a finite positive Borel measure on Q 
and to a signed Borel measure on dQ . 

Proof. Suppose K C VL is a closed set. Cover K with open sets U X (N) 
contained in balls of radius less than 1/N centered at points x G K. Since 
dQ is totally disconnected we may assume this is a cover by clopen sets. Since 
K is compact, a finite subcover K C {U m (N)} exists, and 

K = n™ =1 [u m u m (N)}. 

It follows that the sigma algebra generated by the clopen sets is the Borel 
sets. 

Lemma H21 shows that A 9 (1q, E) is a premeasure on the algebra of clopen 
sets. This premeasure has [TTJ p. 30] a unique extension to a Borel measure 
on Q. 

For E <Z Vt c this argument shows that —A q (ln,E) provides a positive 
Borel measure on Q c . Pasting these two positive measures together yields 
two Borel measures on dQ, the positive Borel measure 

\A q {l n , E)\ = A q {ln, Entt)-A q {l n , En tt c ), 

and the signed Borel measure 

A q {l n , E) = A g (l n , E n fi) + A ? (l n , ^nO c ). (4.3) 

□ 
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The measures fi(E) = A 9 (1q, E) constructed above are finite Radon mea- 
sures [T7J Thm 7.8 and p. 216]. In particular for every Borel set E C dQ, 



fx(E) 



inf{//(0), E CO, O open}, 



and 



f-i(E) = swp{fi(K), K C E, K compact}. 



With the additional assumption that q(x) = for all x in some neighborhood 
of dQ, Theorem 14.41 shows that the positive measure cr(E) = \A q (lf l ,E)\ is 
strictly positive on all nonempty clopen subsets E of Q. 

The following argument shows that the measures /i will be absolutely 
continuous with respect to the measures a. The first step is the next lemma, 
whose proof is straightforward. 

Lemma 4.7. Suppose K C dQ is compact, O C dQ is open, and K C O. 
Then there is a finite collection of clopen sets E±, . . . ,Ek such that K C 



Proposition 4.8. The measures fi = A 9 (1q) are absolutely continuous with 
respect to any Radon measure a which assigns positive measure to each clopen 
set. 

Proof. Using the regularity of our measures, for any Borel set E with fi(E) > 
there is a compact set and an open set M such that C E C M and 
n{E)/2 < /i(^) < n(Op) < 2fi(E), and similarly for cr(E). By taking the 
intersection M fl O a and the union U K a , the same compact and open 
sets may be used for both measures. By Lemma \4. 71 there is a finite covering 
of K by clopen subsets E\, . . . , Ek of O. Since a assigns positive measure to 
all nonempty clopen sets, it follows that o~(O a ) > 0, and so &(E) > 0. 



5 Operator Theory 

Let C(dQ) denote the Banach space of continuous functions F : dQ — > R 
with ||F|| = max xed g \F(x)\. Let M(dQ) denote the Banach space of signed 
finite Radon measures /i on dQ, with 



UE k C O. 



□ 
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As noted above, the measures f-i(E) = A q (1q,E) are finite Radon measures. 
By the Reisz Representation Theorem p2| pp. 216-217] or [231 p. 82], M(dQ) 
is the dual of C(dQ). 

The Dirichlet to Neumann function A q (F, E) maybe extended to the 
Dirichlet to Neumann map A q : C(dQ) — >■ M(dQ). Suppose F : <9£ — > R is 
a continuous function with finite range, or equivalently is a linear combina- 
tion of characteristic functions of clopen sets, F = Xwt=i a fcln(fc)- With this 
domain, define the Dirichlet to Neumann map to be the operator 

n 

A q F(E) = A q (F, E) = Y, «feA g (ln (fe ), E) (5.1) 

k=l 

Since the definition (13. 6p is linear in F : <9C? — > R, the measure valued map is 
consistently defined for clopen sets fl(fc) and in <9(?. The following theorem 
establishes symmetry of the Dirichlet to Neumann function. 

Theorem 5.1. Suppose £1(1) and £1(2) are nonempty clopen subsets of dQ . 
Then 

A q (l m ,£l(2)) = A q (l m ,£L(l)). (5.2) 

Proof. Beginning with some notation, let <f),ip E A, with O<0<1,O<^< 
1, 

nx> lo, xedg\n(i)i' n ] lo, xedg\£i(2)J- 

For 5 > let 

iV(5) = { x eg\ d(x,dG) < 5}. 
For 5 sufficiently small x G N(8) implies 



rl, d(i,n(l))<i, 1 -xfl, 0(2)) < 

lo, d(x,ac\nm) < <u' w lo 



0, d(x,<9£\ft(l)) < 5,J' rv 7 10, d(x,dg\£l(2)) <5, 



The claimed symmetry of the Dirichlet to Neumann function is elemen- 
tary if Q is finite. For j = 1,2, let u^u\ be the g-harmonic extension of lny)- 
By (13. 5 p and the symmetry of the finite graph Dirichlet to Neumann map, 

A g (l n (i),n(2)) = / u n{1) Lil) = ^ 9 u u n{1) {v) = ^2u n{2 )duu n{1] (5.3) 
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The argument proceeds by reduction of the general case to the finite 
graph case. Let Q n be an increasing sequence of finite subgraphs exhausting 
Q. Assume each Q n includes all edges with a point x where <p\x) ^ or 
ip\x) 7^ 0. Also assume that points in the (relative) boundary of Q n are 
contained in N(l/n), with l/n < 5. Define 

fii(n) = dg n n {x I <j>(x) = i}, n 2 (n) = dg n n {x \ ^(x) = i}. 

Let f n be the g-harmonic function on Q n with boundary values F n (x) = 
(f)(x) for x G 9^ n , and let g n be the g-harmonic function on Q n with boundary 
values G n (x) = ip(x). The summands in 



are positive terms, and letting (pi = 1 — there is a sum of negative terms 





Since < f n < 1, these expressions combine to give 




with C<j, independent of n. 
Integration by parts gives 








n 




SO 




The Cauchy-Schwarz inequality then gives the Lipschitz estimate 



\fn(y) ~ fn{x)\ < ^C<pdist(x,y), 
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the distance computed in Q n . Similar estimates hold for g n . 

Since the sequence {f n } is bounded and equicontinuous, it has a subse- 
quence, still denoted {f n }, converging uniformly to / on finite subgraphs of 
Q. The argument in the last paragraph of the proof of Theorem 12.71 shows / 
is g-harmonic. The function / satisfies the same type of Lipschitz estimate, 
the distance now computed in Q. Since / is uniformly continuous on Q, it ex- 
tends continuously to Q. In particular, f(x) = 1 for x G f2(l), and f(x) = 
for x G dQ \ f2(l), so / = /n(i)- Similar comments apply to {g n }, which has 
a limit function 5^(2) • 

Making use of (15. 3p . 

A q {l m ,Q{2)) = [ f m lAj) = lim / f n ty = HmA g (ln l( „),fi 2 (n)) (5.4) 
Jg n Jg n n 

= limA ff (lna(n),fii(n)) = lim / g n L(j) = A g (ln(2), 

n n I r 

□ 

Theorem 5.2. T/ie Dirichlet to Neumann map A q : C(dQ) — > M(dQ) is 
densely defined, symmetric, and nonnegative. 



Proof. As noted in the discussion preceding Lemma [3721 distinct points x, y G 
<9(? are contained in disjoint clopen sets. The Stone- Weierstrass Theorem 
then shows that the linear combinations of characteristic functions of clopen 
sets are dense in C(dQ). 

Let (/i, F) = f g g F dfi denote the dual pairing of measures and continuous 
functions. Suppose F and G are in the domain of A q , 

m n 

j=i k=i 

Then by Theorem 15.11 

(A q F, G) = J2 a iJ G d A i ln F(j) = ^2 a jPkA q (l nF{j) , Q G (k)) 

j j,k 

= ^a i /3 fc A ? (l nG(fe) ,O i ,(i)) = (A q G,F), 
establishing the symmetry of A q . 
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A limiting argument as used in Theorem 15. II will establish the positivity. 
Work in the finite dimensional subspace of functions F which are linear 
combinations of the characteristic functions of a fixed finite collection of 
clopen sets. Following (j5.4p and (15. 3p . for each pair Q(i), of these clopen 



As discussed in Proposition I3.1[ the corresponding finite graph quadratic 
forms (Aqf n , f n ) are nonnegative, making the limit nonnegative. 
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